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1.0  INTRODUCTION 

The  work  consisted  of  writing,  modifying,  and  running  computer  programs 
in  BASIC  and  Fortran  and  using  these  computer  codes  for  theoretical  modelling  of 
the  NRL  theory  of  low  frequency  relaxation  phenomenon. .  Data  analysis  performed 
Included  work  on  dynamics  of  entangled  polymer  melts,  physical  aging  of  glasses, 
dielectric  relaxation  near  the  glass  transition,  and  transient  luminescence, 
transport  and  photoconductivity  in  ehalcogenide  glasses.  Specifically,  the 
following  physical  problems  were  theoretically  modelled  and  compared  with  experi¬ 
mental  data. 

2.0  POLYMER  MELTS 

Application  ot  the  NRL  model  to  polymer  melts  provides  a  quantitative 
explanation  for  the  viscoelastic  properties  of  entangled  polymer  chains  including 
G(t) ,  G" (w) ,  Gn*Jc#  and  n  as  well  as  their  molecular  weight  dependences.  The  model 
predicts  simultaneously  the  temperature  dependence  of  n  as  given  approximately  by 
the  Vogel  form  and  the  shift  factor  a_  in  the  WLF  form.  Modelling  is  quantitatively 
successful  for  entangled  linear  as  well  as  branched  polymers. 

3.0  PHYSICAL  AGING 

Aging  effects  are  found  to  be  contained  in  an  increase  in  the  NRL  model 
parameter  n  as  aging  proceeds,  resulting  in  a  shift  of  the  entire  spectrum  of 
relaxation  times.  The  class  of  creep  compliance  functions  corresponding  to  the 
NRL  model  relaxation  modules  turns  out  to  naturally  describe  Andrade  creep  behaviour 
and  at  long  times  produces  flow  like  behavior.  The  modelling  allows  quantitative 
comparisons  with  creep  measurements  on  glasses,  amorphous  metals,  and  polymers. 

4.0  NaB-ALUMINA 

Frequency  and  temperature  dependences  of  the  dielectric  response  of  the 
fast  ionic  conductor  NaB-alumlna  is  found  to  be  described  by  the  NRL  model.  In 
particular  unusual  properties  in  the  region  from  50  to  300K ,  in  which  a  transition 
from  a  glass- like  to  a  liquid. like  phase  takes  place  now  be  understood  by 

means  of  this  modelling  of  NaSraluoina  data. 

5.0  TRANSPORT  IN  CHALC0GENIDE  GLASSES 

The  chalcogenlde  glasses  are  the  most  thoroughly  studied  amorphous  semi¬ 
conductor  systems.  The  transient  transport  data,  transient  optical  data,  dynamical 
dielectric  relaxation  data,  and  the  volume  and  enthalpy  recovery  data  of  chalcogenlde 
glasses  are  found  to  conform  to  the  NRL  model  predictions.  All  the  available  data 


1 


has  been  analyzed,  in  particular  the  recent  transient  optical  and  transport  data, 
allowing  for  the  first  time  a  clear  distinction  between  predictions  of  various 
models.  NHL  model  behavior  is  found  in  all  cases  and  transport  by  small  polaron 
hopping  is  confirmed. 

6.0  PUBLICATIONS 

The  attached  paper  resulted  from  work  done  on  the  contract. 
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STRESS  RELAXATION  AND  LOSS  MODULI  MEASUREMENTS  IN  LINEAR 
AND  BRANCHED  POLYMER  SYSTEMS  CORRELATED  WITH  STEADY  STATE 
VISCOSITY:  MOLECULAR  WEIGHT  DEPENDENCES 


by 

v 

K.  L.  Ngai  and  R.  W.  Rendell 
Naval  Research  Laboratory,  Washington,  D.C.  2037S 

Nearly  a  decade  ago,  de  Gennes*  Introduced  the  concept 
that  entangled  polymer  chains  relax  by  reptation.  Doi  and 
Edvards  (DE)1  have  developed  a  Model  relating  the  dynamics  of 
reptatlng  chains  to  viscoelastic  properties.  They  assuaed 
that  in  polymer  Melts  and  solutions  with  chain  entangleaents , 
reptation  would  be  the  dominant  notion.  The  stress  relaxa¬ 
tion  modulus  C(t)  in  the  terminal  region  has  been  given  and 
from  it  the  plateau  modulus,  the  viscosity  0  and  the  recover¬ 
able  compliance  J£  can  be  obtained  together  with  their 
dependences  qn^M.  In  particular  rptf  ,  in  close  agreement 
with  the  tpOr  law  of  large  mol .  wt.  linear  polymers. 
Acceptance  of  the  de  Gennes-DE  model  will  follow  confirmation 
of  reptation  by  microscopic  probes  such  as  neutron  scatter¬ 
ing.  Unfortunately,  the  results  of  a  recent  experiment1 
designed  for  such  a  test  fail  to  support  the  existence  of 
reptation.  Until  reptation  is  confirmed  we  must  be  aware  of 
it  being  at  the  present  tine  an  assumption  only.  The  purpose 
of  this  work  is  to  show  that  the  fundamental  mechanism  (FM) 
behind  a  unified  model  (UM)  of  dynamical  processes  in  con¬ 
densed  matter4  readily  provides  an  explanation  for  the  ipT 
dependence,  as  well  as  viscoelastic  properties  including  G(t), 
G"(u>),  G°  and  J~.  The  reptation  assumption  is  not  required. 
Besides, "the  saSe  FM  predicts  simultaneously  the  temperature 
dependence  of  i)  aa  given  approximately  by  the  Vogel  form,  and 
the  shift  factor  a-  in  the  WLF  form.1  The  reptation  model 
does  not  address  tie  temperature  dependence  problem  at  all. 
Moreover,  the  FM  is  applicable  to  both  entangled  linear  apd 
branched  polymers,  while  the  reptation  model  in  its  present 
form  is  inapplicable  to  branched  polymers. 


So  coherent  description  of  the  FM  and  the  UM  will  be 
attempted  here.  .  For  details  the  reader  is  referred  to  Refs. 
4-9.  It  has  been  shown  that  the  relaxation  rate  t . 1  of  a 
mode  i  is  modified  bv  the  FM  to  be  time  dependent  according 
to  t* 1 exp (-n . y) (»  t)  nl  where  y=0.577,  »  is  a  cut-off  fre¬ 
quency  and  n,  is  the  infrared  divergenceexponent  with  OSn  <t. 
The  rate  equation  for  mode  relaxation  is  modified  from  iQjJdt* 


-(l/YjJQj  to  dQi/dt=-(l/T.)exp(-niY)  (u»ct)  i  Qj.  The  solu¬ 
tion  of  which  is  Q.ttJsQ^OjexpI-expC-njYJt1  ni/(l-ni)liu»cI1il3 

Q. (0)expl-t/t  .)*  i,  which  for  n-»0  reduces  to  Q.(0)exp(-t/t.). 
Wnen  n#0  it  is  a  fractional  exponential  and  is  shifted  to1 


rpi  =  lexp(niY)(l-ni)tiwcni),/(1'ni)  (1) 


This  specific  .fractional  exponential  relaxation  function 
being  a  consequence  of  a  FM  has  led  us  to  work  out  its  unify¬ 
ing  predictions  in  many  fields  and  Mterials.  We  point  out 
here  the  aforementioned  shift  front,  to  t  .,  which  we  shall 
in  this  work  use  to  derive  the  rp4T  relation,  etc.,  has  been 
repeatedly  tested  and  verified  in  other  areas  including  elec¬ 
tron  transport4*8  in  amorphous  SiO.  and  chalcogenide  glasses, 
ion  transport  in  conductors,*  dielectric  and  structural  relaxa¬ 
tions  near  7  lor  polymers  and  smII  molecule  glasses,7  y 
relaxation  of*bulk  polycarbonate,8  volume  recovery*  of  PVAc, 
etc. 


To  start  with  consider  the  Rouse  sMdel.  The  Rouse 
internal  nodes  are  periodic  functions  of  bead  positions  along 
the  chain.  The  1-th  mode  has  characteristic  period  equal  to 

N/i  and  relaxation  tine  /(iriVkT  Mr  where  If  is 

1  o  o  o 

the  iMitoner  mol.  wt.  and  t  is  the  friction  coefficient.  The 
lower  nodes  govern  the  largo  scale  motions  and  have  the 
longest  relaxation  tines.  The  Rouse  model  gives  a  good 
description  of  data  on  linear  chain  polysters  of  low  mol.  wt., 
with  rp4t  for  H<H  .  For  longer  chains,  i|  is  a  much  stronger 
function  of  H.  cFor  M>M  ,  the  rpdr*  power  dependence  holds. 
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In  some  linear  system*  |i  appears  to  be  larger  than  3.4.  In 
branched  systems,  |i  can  be  considerably  larger  especially  at 
high  mol.  wts.  The  change  in  M  dependence  occurs  smoothly  as 
M  varies.  When  H  exceeds  H  a  plateau  region  appears  in 
J(t),  C(t),  C'  (u»)  and  C"(u»)?  These  H-dependences  are  a 
direct  indication  of  chain  entanglement  which  sets  in  when 
M>H  and  gives  rise  to  a  temporary  network  structure  due  to 
intirlooping  of  the  chains.  Thus  the  onset  of  entanglement 
introduces  a  coupling  between  chains.  If  H  is  the  average 
mol.  wt.  between  entangelements  commonly  defined  through  the 
plateau  modulus  by  H  «pRT/G°,  the  number  of  enganglement 
couplings  per  molecuf e  is  M/M -1.  Internal  Rouse  nodes  with 
i>N/2M  Si  involve  cooperative  notions  over  distances  smaller 
than  tfle  entanglement  spacing.  They  are  of  sufficiently 
short  range  and  are  unaffected  by  entanglement  coupling.  In 
our  model  this  implies  the  coupling  of  the  i>i  Rouse  nodes 
to  the  correlated  states4  •  in  the  polymer  melt  or  solution 
cannot  occur  through  the  entanglement  junctions  and  it  is  so 
weak  that  effectively  a.=0  for  all  i>i  and  they  are 
unshifted.  On  the  other  hand,  lower  mode  smtions  with  i<i  , 
which  impose  cooperations  that  extend  beyond  the  entangleamnt 
spacing,  are  coupled  via  the  entanglement  junctions  to  the 
correlated  states  of  the  polymer  system  and  have  nonzero  n. 
values.  These  nodes  are  slowed  down  and  according  to  Eq.  (1) 
their  T.'s  are  shifted  to  the  large  t  .  values.  This_is 
because1!!!  x.»l  which  follows  from  the  expected  size4  *  of  w 
and  of  the  typical  relaxation  tines  in  the  plateau  region  anS 
the  terminal  zone.  The  degree  of  cooperation  ia  largest  in 
the  motion  of  the  terminal  relaxation  mode  i=l  and  decreases 
ss  i  increases.  Hence  n.  is  maximum  for  i=l  and  decreases  as 
i  increases  towards  i  .  1The  "gap"  in  the  relaxation  time 
spectrum  between  the  other  i<i  modes  and  the  i=l  mode 
creates  a  plateau  region.  The  gap  is  not  a  true  gap,  having 
some  contributions  from  the  higher  i  modes  in  the  group  i<i 
that  have  smaller  n.  and  thus  smaller  shifts  from  i.  to  x  .. 

In  other  words,  since  n.>a.>  ...  >n.  -1,  the  shift  tlog  p 
x  .-log  X.)  is  not  uniform  with  respect  to  i.  The  transition 
rEgicn  comprises  an  unshifted  Rouse  spectrum  of  i>i  modes. 
Each  of  the  Ki  modes  is  individually  shifted  and  modified 
to  have  the  fractional  exponential  form  of  exp(-t/x  .)  i. 

The  X  .  i.ow  characterizes  the  time  scale  of  the  i-tn  non- 
exponl&tial  relaxation  mode.  The  spacings  between  the  T.'s 
are  different,  much  more  nonuniform  and  wider  than  the  p 
rigidly  shifted  Rouse  spectrum.10  This  can  be  seen  from 
Eq.  (I)  as  due  to  the  (1-n.)  *  power  and  that  n  >n«>  ••• 

>n  -1.  The  resulting  spectrum,  even  though  individual  com¬ 
ponents  are  nonexponential,  is  much  narrower  than  the 
rigidly  shifted  Rouse  spectrum  as  we  shall  see  later.  Recall 
that  T.«dl  for  all  i.  However,  for  i<i  ,  T«  is  shifted  to 
X  .  which  now  varies  as  lr'v  "i  .  Since  2/(l-n^)R2  ve  see 
tie  pseudo-gap  or  the  plateau  region  widens  as  (^increases. 

The  separation  between  I  ,  and  t  ,  is  much  larger  than  that 
between  x.  and  X,.  Hence ‘the  terminal  region  is .almost 
entirely  due  to  f  ..  X  .  is  proportional  to  (T/l  K«trl. 

Its  contribution  to  viscosity  q  according  to  the  Rouse  model 
is  (pRT/M)Tpl«w,l  .  That  the  terminal  relaxation  time  X^j 

and  q  have  different  H-dependences  is  not  in  sgreeswnt  with 
experiment.  This  last  failure  is  not  surprising  because  the 
Rouse  model  should  not  be  applicable  in  its  original  form  at 
high  mol.  wts.  Entenglesttnt  means  the  chain  cannot  execute 
all  Rouse-like  motions.  Some  special  confinement  such  as  the 
-  fixed  tube  used  by  DR  ia  appropriate  though  the  total  suppres¬ 
sion  of  sustained  transverse  motions  ia  the  tube  model  is 
unrealistic.11  Nevertheless  we  consider  next  the  chain  ia  the 
tube  sad  rederive  all  the  results  with  eur  IN.  The  reader 
should  realise  that  the  tube  model  referred  to  here  means 
only  chain  confinement  by  a  tube  and  no  captation  is  invoked. 
In  fact  the  thrust  of  this  work  is  to  shew  that,  without  the 
reptatioa  hypothesis,  the  FH  which  has  already  explained  many 
properties  in  other  related  and  unrelated  fields,  once  mere 
satisfactorily  accounts  for  the  polymer  melt  properties  when 
the  chain  is  confined  in  a  tube. 

The  development  initially  follows  closely  the  works  of 
DE  and  of  Graessley11  except  we  replace  the  reptatioa  modes 
by  our  FH  shifted  fractional  exponential  relaxing  Rouse  modes 
with  shifted  t.  for  Mi  .  Configuration  rearrangements  in 


t 
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the  terminal  zone  are  accomplished  essentially  by  the  1=1 
mode.  The  relaxation  modulus  is  given  by  G(t)=C^exp- 
(t/t  j)'  “l  where  Gj.  is  the  plateau  modulus  given  by  4vNfcT/S 
and  independent  of  H,  The  viscosity  n  and  related  to  G(t) 

by  n  =  J^C(t)  dt  and  J°  =  r|*^J^tG(t)dt  can  be  calculated.  We 


obtain  r) 


"  C».(1‘ 


ni> 


r(l/(l-n,))Tpi  and 


(l/Gpd-nj) 


r(2/(l-nI))/r2(l/(l-n1)),  where  T  is  the  Ganaia  function.  For 
n=0.5,  the  product  J4c£=3,  t  ^  and  0  both  «M*.  For  n=0.4, 

J®0?=2.04,  x  .  and  0  both  The  transition  region  again 

consists  of  fie  unshifted  (i.e.  n  =0)  Rouse  nodes  for  i>i  . 
The  plateau  region  has  some  contribution  from  the  unevenly 

shifted  (i.e.  m.>a.+.)  Rouse  modes  i=2 ,  3  .  i  -1  in  a 

tube.  The  derivation  of  a  rheological  constitutive  equation 
within  the  tube  model  by  DE,  and  the  predictions  of  stress- 
strain  experiments  which  follow  from  this,  is  largely  unaf¬ 
fected  by  our  replacement  of  reptation  dynamics  by  FM  dynam¬ 
ics.  This  is  due  to  the  approximate  factorizability  of  the 
time  dependence  and  the  strain  dependence.  Thus  many  of  the 
successful  rheological  predictions  of  OE  are  carried  over 
into  our  model  while  we  also  retain  the  advantages  of  the  FM 
dynamics.  We  present  finally  the  predictions  of  our  model 
for  both  linear  snd  branched  polymers.  The  results  are  dis¬ 
played  in  Figs.  1-4.  The  G"(ui)  data  and  the  predictions  are 
both  for  the  terminal  zone.  Note  that  for  narrow  distribu¬ 
tion  linear  polystyrene  melt  (Fig.  1),  the  G"(u>)  dispersion 
is  well  fitted  by  model  prediction  with  n=0.42S  which  at  the 
same  tine  predicts  both  t  .  and  q  should  be  <4T  with  p  = 
2/(1-0.425)23.5  in  agreement  with  the  observed  M3'4  behavior. 
For  branched  polymers  (Fig.  2),  the  fit  to  G"(v>) -demands  a 
larger  n,  value  of  0.7  which  predicts  that  n“M°' '  again  in 
agreement  with  data  in  these  branched  polystyrenes  at  high 
nol.  wts.  The  linear  relaxation  modulus  G(t)  curves  are  our 
model  predictions  obtained  from  an  unshifted  Rouse  spectrum 
for  i>i  ,  and  for  iSi  we  chose  n.  to  rise  linearly  from  zero 
value  at  i=i  to  n.  as  i  decreases  towards  unity.  The  fit  to 
the  G(t)  data  of  Fig.  3  for  concentrated  solution  of  poly¬ 
styrene  in  Aroclor  1248  requires  n.  =  0.475  which  implies 
that  ryUT  with  p  =  2/(l-0 .475)=3.8.  On  the  other  hard 
branched  polystyrenes  (Fig.  4)  of  lower  mol.  vt.  than  that  in 
Fig.  2  requires  n.  ~  0.56  to  0.61  which  predicts  p  ~  4.5  to 
5.1.  Our  model  predicts  that  dispersion  of  the  viscoelastic 
functions  (determined  by  n.)  at  the  terminal  region  is  cor¬ 
related  with  the  ncdr  behavio 


ior  via  the  relation  p=2/(l-n1> 


and  the  size  of  J°G°  via  the  relation  J°G°=(l-n  )r(2/(l-n, ))/ 
a  c  n  C  A  1  1 

r^(l/(l-n.)).  For  linear  polymers,  typically  we  find  from 
dispersion  that  nj  ranges  between  0.4  to  0.5  for  H>Mc>  These 

n.  values  give  3.3SpS4  and  2SJ°G°S3,  all  in  agreement  with 
data.  For  branched  polymers,  as "determined  from  dispersion, 
n.  typically  has  larger  values  which  implies  larger  p  and 

J  cj.  When  n,  is  a  function  of  M,  the  imdel  further  predicts 
that  both  0  and  become  M  dependent.  All  these  features 
are  in  remarkable  agreement  with  data.  Details  and  further 
implications  of  the  model  will  be  presented  elsewhere. 

Sachs/ Freeman  Assoc.,  Inc. 
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